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Abstract 

Air transport is a key infrastructure of modern societies. In this 
paper we review some recent approaches to air transport, which make 
extensive use of theory of complex networks. We discuss possible net- 
works that can be defined for the air transport and we focus our atten- 
tion to networks of airports connected by flights. We review several 
papers investigating the topology of these networks and their dynamics 
for time scales ranging from years to intraday intervals, and consider 
also the resilience properties of air networks to extreme events. Finally 
we discuss the results of some recent papers investigating the dynam- 
ics on air transport network, with emphasis on passengers traveling in 
the network and epidemic spreading mediated by air transport. 
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1 Introduction 



Transport systems, along with other infrastructures hke power grids 
or communication networks, are fundamental elements of our soci- 
eties and economies. They guarantee the high level of mobility that 
we all experience, and which is vital for the cohesion of markets and 
for the quality of life of citizens. Moreover, transport systems enable 
socio-economic growth and job creation. When such fundamental in- 
frastructures experience a random failure or are intentionally targeted 
by terrorist attacks, the whole society is severely affected. The air 
transport system is no exception. In 2008 this industry generated 32 
million jobs worldwide, of which 5.5 of them were direct, and con- 
tributed with USD 408 billion to the global gross product [T]. On the 
other hand, its vulnerability and the consequences for citizens' mobil- 
ity clearly appear when a strike or the eruption of a volcano interrupt 
the normal behavior of the system [21 [3]. 

The theory and application of complex networks has experienced 
a tremendous growth in the last decade [H [5] . In spite of its young 
age, the great variety of tools developed for the analysis of different 
topologies 1^ has favored a better understanding of the structure and 
dynamics of many real- world systems [7]. Remarkably, the complex 
networks approach has explained the appearance of emergent phenom- 
ena in many systems composed by a large set of interacting elements. 
Well known examples include social systems, e.g. the study of net- 
works of acquaintances or the diffusion of contagious diseases [8] , the 
Internet [9j , and applications to neural dynamics [101 [TT] . It is not sur- 
prising that the complex network methodology has been successfully 
applied to different transportation modes, including streets [T2| [T3]. 
railways [H] , or subways [HI [16] . 

In this paper we present a review of the literature related to the 
application of complex network theory to the air transport system. As 
it will be clear, several problems have been investigated so far. For 
instance, the description of the topological and metric structure of the 
network is of great importance for understanding the business strate- 
gies adopted by different airlines, for assessing passengers' mobility in 
the presence of direct and indirect connections, or for investigating 
the time evolution of air transport, while it adapts to changes in the 
passengers' demand and reacts to economical external forces, such as 
deregulation. Another aspect of interest is the dynamics taking place 
on the network. A paradigmatic example, reviewed in this paper, is 
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the spreading of infectious diseases worldwide and the role that air 
transport has in enhancing the speed of epidemic propagation. 

The future presents many challenges for the air transport system 
and complex network theory is likely to play a more and more sig- 
nificant role in tackling these challenges. First of all, air transport 
is increasing worldwide at a very fast pace. Policy makers are aware 
of the fact that the current system will be at its capacity limits in 
few years because of the increase of traffic demand and new business 
challenges. For this reason, large investment programs like SESAR 
in Europe and SingleSky in the US have been launched. Also, policy 
makers have stressed the importance of fostering the resilience of the 
system, and of its capacity of recovering the required mobility after 
an external shock [17J. Moreover the future will require an increas- 
ing degree of integration among different transportation modes. This 
problem finds a natural description in terms of multi-layer represen- 
tation of complex networks [H]. Clearly, all these issues are not only 
relevant for the air transport itself, but they have important implica- 
tions for the society whole. 

This review is organized as follows. Section [2] describes the main 
components involved in the air transport system, which are the basis 
for the construction of different network representations. Section [3] re- 
ports the most important facts about the topologies of such networks, 
including their dynamical evolution, and the models that have been 
developed to explain these characteristics. Section U] reviews the main 
dynamics that have been studied on top of this network and Section 
[5] discusses the role of network topology in the resilience and vulner- 
ability of air transport system. Finally, Section [6] draws some final 
conclusions, and presents some open lines of research. 

2 Networks for the air transport 

Many complex systems can naturally be represented by one or more 
networks. For instance, the Internet can be represented as a set of 
nodes (the webpages) connected by links (the hyperlinks) [19l [20] . 
The same system can also be represented by considering the routers 
as nodes and their physical connections as links |21] . 

This multiple network representation property is shared also by 
the air transport system, and therefore one should first decide which 
network is investigated. The air transport system is composed of a 
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large number of different elements, interacting and working together. 
The mobility of passengers is just the final result and it is clearly 
of high importance from a social point of view. Therefore it is not 
surprising that most analyses have been focused on the mobility of 
people, disregarding other technical details. When this point of view 
is followed, the construction of network is straightforward. Nodes 
represent airports and a link between two nodes is created whenever 
there exists a direct flight between the two airports associated with the 
nodes. From this point of view the airport network is the projection of 
the bipartite network, whose first set of nodes is composed of airports, 
the second set of flights, and a link exists between a flight and an 
airport if that flight departs or arrive to that airport. Clearly in 
this way additional sources of information, like scheduling, types of 
flights, or airlines, are disregarded. The projected network of airports 
is naturally a directed graph, where, in general, two directed links can 
exist between two nodes A and B, one describing the flights from A 
to B and one from B to A. The projected network has also a natural 
weighting scheme, given by the number of flights that are present 
(in the investigated time period) between the two airports. These 
networks are termed flight networks. 

Additionally, not all flights are equivalent. The number of available 
seats in each aircraft can be very different, ranging from the 50 pas- 
sengers of a small regional jet, up to the 853 seats of an Airbus A380. 
As a consequence, it has been proposed the association of weights to 
links, proportional either to the frequency of connections or to the 
number of transported persons - see Section 13.21 As usual, from a 
weighted directed graph one can construct other graphs by neglecting 
information. For example, taking the difference between the weights 
from ^ to and from B to A, one can construct a directed network 
where only one directed link exists between two nodes. By neglect- 
ing the weights one can obtain an unweighted network, where only 
topology matters, and by neglecting also directionality one can obtain 
simple binary graph. All these alternatives have been investigated in 
the literature. In any form, the flight network is probably the most 
investigated network in air traffic studies. Here we will review mostly 
the properties of these networks in the following of this paper. 

However while complex networks are traditionally considered as 
static objects, it is clear that time is an important feature of any 
movement. This is especially true in the case of the air transport, 
because passengers may have to use several flights to get to their 
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destination. If one considers only static representations of the network, 
there is no way of knowing the real dynamics of passengers, i.e., if 
one needs to wait 2 or 10 hours in a airport before taking the next 
connecting flight. In Section H] we will review some solutions that have 
been proposed to investigate the problem of indirect connectivity of 
passengers. 

Another important aspect of flight networks is that they can be 
naturally decomposed into many subnetworks. For example, flight 
networks can be decomposed by considering separately one subnet- 
work for each airline. The analysis of the networks corresponding to 
a single airline has been performed, for instance, in Refs. |22| 123] . 
To the best of our knowledge, no paper has been published on the 
analysis of interdependencies between the subnetworks corresponding 
to different airlines, or of different alliances of airlines. We expect 
this topic will gather an increasing interest, especially because it can 
be considered a special case of a multi-layer representation of com- 
plex networks [18j. This framework would also allow the study of the 
relations between the air and other transportation modes [241 125j. 

It is important to stress that networks different from the flight 
networks can be constructed and are likely to be very important for 
the understanding and modeling of air transport. The structure of the 
airspace is one of these cases. Nowadays aircraft do not travel along 
the straight line (geodesic) connecting the departure and destination 
airport. On the contrary, they must follow some fixed airways, defined 
as union of consecutive segments between pairs of navigation aids 
(navaids in short). While such constraints are actually imposed in 
order to improve safety (as it is easier to control ordered flows of 
aircraft) and capacity (the workload of controllers is reduced, and 
thus they can control a higher density of aircraft), bottlenecks may 
appear in some central zones of the airspace, or where several busy 
airways converge. Navaids can be used to create a network, where 
nodes are navaids, and links represent airways. To the best of our 
knowledge, only very few studies have considered this type of network. 
For example, Cai and coworkers [26] investigated the Chinese air route 
networks. 

Another example of networks considers reactionary delays and their 
effect on passengers. These are situations in which a flight cannot take 
off on time because of a delay in another flight. This occurs, for in- 
stance, because of the late arrival of the aircraft, or of the crew itself 
[27]. Different networks may be created to study this phenomenon. 
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Nodes may represent crews, with a link between them when they share 
the same aircraft. Alternatively nodes may also represent airports, 
connected whenever the same aircraft has to serve them in a sequen- 
tial fashion. The identification of the central elements of these graphs 
may help in highlighting the critical points for the dynamics of the 
system, and would thus allow the creation of better mitigation strate- 
gies. To the best of our knowledge, this topic has not yet been studied 
within this framework. 

Finally we mention a very recent approach to construct networks 
of air traffic safety events [29]. When two aircraft are too close, an au- 
tomatic alarm, termed Short Term Conflict Alert (STCA), is activated 
and the air traffic controller is supposed to give instructions to the two 
pilots in order to avoid a collision. One important question is whether 
STCA are isolated events or whether aircraft initially involved in a 
STCA are likely to be involved in other STCAs with other aircraft in 
the near future and so on, creating a cascade of events. This possibil- 
ity signals the fact that the controller suggests a local solution without 
forecasting unintended consequences of her instructions. By using a 
dataset of automatically recorded STCA, authors of [29] mapped this 
problem into a network of STCAs which in turn can be mapped in a 
network of aircraft, where two nodes (aircraft) are connected if they 
were involved together in a STCA. These networks shows topological 
regularities and might shed lights on the aircraft conflict resolution 
dynamics. 

3 Topological analysis 

3.1 Unweighted air transport networks 

The analysis of the structure of the flight network in air transport, 
especially when focused on individual airlines, began years before the 
formalization of the complex network theory. This type of analysis 
was motivated by the aim of deflning the most efficient structures of 
flights for a given airline \28\ I30j . both in terms of yields (and, thus, 
proflt) and of passengers' mobility. The proposed solutions can be 
grouped into two classes: 

Point to point: in this configuration, a different aircraft serves each 
pair of airports in the network, or at least those pairs where the 
passengers' demand is enough to justify the connection. While it 
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has the advantage of offering direct connections to all passengers, 
it also requires a high number of aircraft for covering all the pos- 
sible routes. For a completely connected network the number of 
connections increases with the square of the number of airports. 
This strategy of connections was common in the United States 
before the '70 deregulation [31j, and nowadays is still used by 
several low-cost airlines - see left panel of Fig. 1. 

Huh- and- spoke: in this case, connections are structured like a chariot 
wheel (or a collection of such structures), in which all traffic 
moves along spokes connected to the huh(s) at the centre. While 
most passengers must take (at least) two different flights to reach 
their destination, this strategy presents several benefits for the 
airline. In fact, a lower number of aircraft is required, flights 
usually have a higher occupation rate, and the expansion of the 
network to a new airport only requires one new additional flight 
|32ll33j . Today the hub-and-spokes configuration is used by most 
major airlines all around the world - see right panel of Fig. 1. 

While earlier studies were mostly theoretical, the possibilities of- 
fered by the analysis of real systems through the complex networks 
methodology, and the ever-increasing computational capabilities of 
modern computers, have enabled a better understanding of the struc- 
ture of real air transport networks. It is interesting to notice that 
some network characteristics have been confirmed in all studied net- 
works. One important aspect of flight networks is the the fact that 
they show the scale-free feature. This implies the presence of few hubs 
with a very high number of connections, confirming the predominance 
of a hub-and-spoke topology. An example can be seen in Fig. 2, 
which shows the cumulative probability distribution of degree of all 
European airports, considering only internal flights (i.e., flights whose 
origin or destination is outside Europe are disregarded). The right 
panel of Fig. 2 shows a zoom of the extreme tail of the distribution 
and it is clear how few airports have direct connections with a large 
number the destinations in the network, performing a hub function. 

In the literature we surveyed, different methodologies have been 
used for creating and analyzing flight networks. In Tabled] we report 
the values of classic complex network metrics, for the air transport 
system of different countries and considering an unweighted represen- 
tation. The meaning of these metrics is reported below: 

7: in scale-free networks [45] the asymptotic behavior of the node de- 
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Table 1: Example of different topological metrics of flight network, as re- 
ported in several research papers. The asterisk in the Links column indicates 
that the number refers to the number of flights, while in all the other cases 
the column reports the number of connections. 
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Figure 1: Representation of the networks corresponding to Ryanair (Left) 
and Lufthansa (Right), as at 1^* of June 2011. Notice how the network of 
Lufthansa, which is a major airline, is centered around few main airports 
(hubs), while the structure of Ryanair, the biggest European low-cost com- 
pany, has a densely connected core. 
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Figure 2: Cumulative probability distribution of degrees of European air- 
ports. The network has been constructed considering only commercial (both 
regular and charter) flights operated between two European airports the 1^* 
of June 2011. (Right) Zoom of the right part of the distribution. The hubs 
of the network are indicated by red circles. 



gree distribution has a functional form P(k > x) ~ x"'*'. It has 
been pointed out that the real degree distribution of the world- 
wide flight network is a truncated power-law, i.e., it is asymptot- 
ically better explained by the function P{k > x) oc f{x/km), 
where / is 3>n Gxponenticil trunccition function and. h^n is a trunca- 
tion parameter. The values reported in Refs. [Ml EH] correspond 
to the exponent 7. 

7b: the betweenness of a node is a centrality measure quantifying how 
important is a node for movements inside the network. Node be- 
tweenness is defined as the proportion of shortest paths, among 
all possible origins and destinations, that pass through a node 
|46j . The exponent 7s is the exponent of a power law fit of 
betweenness distribution. When the distribution of centrality is 
asymptotically a power-law function, a high exponent 75 indi- 
cates that few nodes are responsible for the efficient routing in 
the network. 

L and L^and' L is the mean length of shortest paths between pairs 
of nodes of the network, i.e., 
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where i and j are two nodes of the network, n the number of 
nodes, and dij the length of the shortest (topological) path be- 
tween nodes i and j. The value of L is usually compared with 
Lrandi that is the mean value obtained in different networks that 
have the same number of nodes and links, but a completely ran- 
dom structure. These random networks are also called Erdos- 
Renyi graphs [1^ . It is worth noticing that the Table shows how 
L is often lower than the corresponding Lrand: indicating, as ex- 
pected, that air transport networks are engineered to efficiently 
reduce the number of connections needed by passengers. 

C and Crand- the clustering coefficient C, and its randomized coun- 
terpart Crandj measures the number of triangles that can be 
found in the network It assesses the probability that two 
nodes, which are connected to a third node, also share a direct 
connection. Similarly to Lj-and, Grand corresponds to the mean 
clustering coefficient of an ensemble of randomized networks. 

The reader may easily notice how obtained values are very het- 
erogeneous. For instance, the exponent of the degree distribution 7 
varies from 1.0 up to 4.161, and the clustering coefficient C from 0.07 
to 0.738. This variability is mainly due to two factors. Firstly, there 
are important differences in the method used in the construction of 
the networks, which are usually not fully explained in the papers. 
The time window represented by the network may be not reported 
[22I [23l m] , and no details are given about the types of flights con- 
sidered (regular passengers flights, charters, cargo flights). Secondly, 
most of the researches have investigated national networks, covering 
few tens of airports. It is well known that some complex networks 
properties, such as, for instance, the scale-free distribution of degrees, 
are meaningful and can be correctly assessed only in large networks, 
where finite size effects are negligible |48j . Moreover the degree of net- 
work heterogeneity is very different if one considers a regional airport 
network or the worldwide network. 

3.2 Weighted network analysis 

As explained in Section [3. 11 the analyses described above are based on 
unweighted projections of the air transport system, that is only the 
existence of direct connections between pairs of nodes is taken into 
account. On the other hand, it can be expected that the structure 
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of frequencies of flights may unveil interesting information, especially 
related with the main routes of movements chosen by passengers. 

Table [2] shows the values of some metrics obtained for different 
weighted networks [52J. When a link between two nodes i and j has a 
weight Wij, it is possible to calculate a weighted version of the degree 
of a node, called strength, as Si = Yl,j'^ij- Notice that the variability 
observed in the metrics of the unweighted networks is amplified in 
the weighted network, because several variables can be used to define 
the value of Wij. For example one can consider the number of flights, 
the number of offered seats, or the number of passengers transported, 
obtaining different weighted networks. 

The definition of the metrics shown of Table [2] is here reported: 

j3: the relation between the strength s (number of flights) of each 
node and its degree k (number of connections) is typically well 
fitted by a power law s{k) ~ A;^. This relation unveils relevant 
information about how capacities are distributed through the 
airport network. 

if one is interested in the assessment of the centrality of airports 
from the point of view of passengers' movements, it is possible to 
relate the strength of a node with its betweenness, i.e., s{h) ~ b^^ . 

9: in order to check whether there is a relation between the frequency 
of connections between two airports, and their connectivity (in 
terms of number of destinations directly served), the weight of 
each link has been related with the degrees of connected nodes, 
leading to Wij « {kikjY . 

The main conclusion that can be drawn from Table [2] is that there 
exists a strong correlation between the degree of a node, and the quan- 
tity of flights and passengers going through it. This fact is in agree- 
ment with the hub-and-spoke structure of the network. The more 
connections a node has, the more passengers will use that node to 
reach other destinations, and thus the frequencies of such connections 
strongly increase. 

3.3 Short and long term evolution of the net- 
work 

Although the scheduling of flights of regular airlines should be defined 
and published almost one year before the actual day of operation. 
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Table 2: Topological properties of different weighted flight networks. 



Country 


Weigth 
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Worldwide 


Available seats 


1.5 


0.8 0.5 


[3511361119] 


US 


Number of passengers 


1.8 




[38] 


India 


Number of flights 


1.43 




m 


4 European airlines 


Number of flights 


(1.06-1.18) 




[50J 


China 


Number of flights 




— 0.5 


im 


Europe 


Number of flights 


1.39 




[51] 



the air transport system is far from being a static structure. On the 
contrary, airlines constantly adapt their scheduling to changes in the 
passengers' demand, both on the short and long term. 

On the short term, the network evolves to answer to the different 
needs of two groups of passengers, namely those moving for work 
and usually traveling from Monday to Friday, and those moving for 
leisure, traveling mainly during the week-end. It is well known that 
the first group is less sensitive to price, but highly values short travel 
durations. The result is that, during the week-end, the number of 
flights is reduced, and the network assumes a more star-like shape 
around the airports of cities of tourist interest. Fig. 3 reports the 
evolution of two metrics, namely the mean degree (k) and the exponent 
of the power-law fit of the degree distribution 7, for the Austrian and 
Chinese air networks at different days of the week. 

A long term adaptation of the network has been also studied in 
several papers. While short term variations are easy to predict, long 
term oscillations in the demand are more complicated to forecast, be- 
cause they are the result of significant changes in macro-economical 
factors, of the competition between different airlines, and of the limi- 
tation imposed by policy-makers. 

The evolution of the structure of the network due to changes in 
the regulations has been extensively studied in the case of the Eu- 
ropean deregulation started in year 1986 [53] (see also Section 1^]) . 
While several papers have focused on the aeronautical and economical 
consequences of the deregulation [Ml Ell [56] , only a few have applied 
complex networks concepts to this problem. Specifically, in Ref. [57] . 
the evolution of the European air network is studied between years 
1990 and 1998, with one network created for each year. The dynam- 
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Figure 3: Evolution of the mean degree (k) (Left) and of the exponent of 
the degree exponential fit 7 (Right) for the Austrian (black squares and solid 
lines) and Chinese (blue triangles and dashed lines) air transport networks, 
through different days of the week. Information about the Austrian and 
Chinese networks have been obtained in Refs. [22] and [23| respectively. 



ics of some basic network metrics have been investigated. The most 
important of them are the classification of airports according to their 
number of connections and the evolution of the number of small and 
big (hubs) airports through time. Results identify two different dy- 
namical patterns. On one side, medium size airports have attracted 
most of the intra-European traffic, creating specialized internal hubs. 
On the other side, the intercontinental traffic was also concentrated, 
but on different hubs, usually big airports. In summary, the global 
structure of the network has promoted an hub-and-spoke structure, 
but with different hubs for internal and intercontinental flights. 

An interesting case study is represented by China. Since the eco- 
nomic reforms started in 1978, which slowly reduced the control of the 
State over the economy, the size of the air network and the number 
of passengers transported have increased at fast pace [58] (see Fig. 
4). Specifically, the number of passengers grew at an average annual 
rate of 17%, mainly driven by business and tourist motivations. As a 
consequence of this expansion, the number of airports also increased 
from 69 in 1980 to 137 in 1998. As forecasted by the classical theory 
[3T] . the structure of the network changed from a point-to-point to a 
hub-and-spoke configuration. Nevertheless this hub-and-spoke topol- 
ogy presents some interesting specificities [59]. First of all, there are 
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Figure 4: Evolution of the number of airports in the China's air network, and 
of the number of passengers from 1980 to 2005. Data correspond to Table 1 
ofRef. m. 



three main hubs, corresponding to the headquarters of the three main 
national airlines, namely Beijing (Air China), Shanghai (China East) 
and Guangzhou (China South). The relative importance of these air- 
ports has significantly changed over time, reflecting the changes in the 
international scene of the corresponding cities. While Beijing was the 
main hub in 1990, it passed its role to Shanghai in 2005 [60] . Finally, 
these three airports form a peculiar triangular subnetwork, accounting 
for a 37.3% of the transported passengers. These characteristics are 
mainly the result of regional economical and social inequalities within 
the Chinese territory, and of China's involvement in the world econ- 
omy. After year 2000, the main topological properties of the network 
have remained stationary. However several changes have occurred, 
mainly related with the appearance and extinction of small airports 
and routes between them |61j . 

To conclude this Section on the evolution of air networks, it is 
worth noticing Ref. [62], in which the evolution of the Brazilian air 
network between years 1995 and 2006 is studied. Similarly to the case 
of China, Brazil experienced from an important growth in its economy. 
Within this period, the total number of passengers transported in 
Brazil increased from 18 millions to 43 millions. Nevertheless, and 
in opposition to what found in the case of the Chinese network, the 
number of airports served by major airlines strongly decreased (211 to 
142), along with the mean degree (13.19 to 10.28). Therefore, while 
the traffic increased, the network reduced the number of connections 
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between secondary airports, getting closer to a pure hub-and-spoke 
configuration centered around Sao Paulo, Brasilia and Salvador. 

4 Dynamics on the air network 

Up to now we have considered the topological and metric properties 
of the flight network, by discussing also the dynamics of the network 
over short and long time scales. However the air traffic network intrin- 
sically describes the space where something, for example passengers 
or goods, moves. Therefore it is natural to consider the problem of 
characterization, modeling, and control of the dynamics on the flight 
network. In this Section we review some recent approaches to dynam- 
ics on air transport network, focusing our attention to the problem of 
the dynamics of passengers when connection is not direct, the emer- 
gence of traffic jams, and the propagation of epidemics through the 
air traffic network. 

4.1 Indirect connectivity and passengers dy- 
namics 

Probably the most important entity moving on the air transport net- 
work are passengers. The structure of the air network strongly affects 
the capability for a passenger to reach destination B starting from des- 
tination A in the shortest possible time and in the most direct way. 
The transition from the point to point system to the hub and spoke 
system had several effects on passengers [32l[33]. On one side, it has 
been argued that the hub and spoke system increased airline's effi- 
ciency and therefore, in a competitive market, it lead to lower prices. 
On the other side, it is clear that spoke cities risk being marginalized, 
being connected in a more indirect way to the rest of the system. The 
emergence of low cost carriers, thanks to the deregulation of the mar- 
ket, created cheaper opportunity of flying and created additional hub 
and spoke subnetworks, where hubs are different from the one used by 
main airlines |63j . These qualitative considerations on the effect of the 
hub and spoke system on passengers have been investigated recently 
in a more quantitative way by making use of network theory, as we 
will show in the following. 

Many papers have focused on the dynamics of the air transport 
network toward an higher spatial concentration, which means that the 
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topology of the network becomes more and more similar to a collection 
of star like structures. For example, several authors have investigated 
the change of topology of the American [M] and European |571 |65] 
air transport network during the transition from the point to point to 
the hub and spoke structure. However as pointed out, for example, in 
Ref. |66) a hub and spoke network requires a concentration of traffic 
both in space and in time. Temporal concentration means that flight 
schedule must be organized in such a way to allow passengers to travel 
between two (or more) spoke cities in a relatively short time, avoiding 
thus long waiting time at the hubs. 

In Ref. j66] authors considered the airline network configurations 
in Europe between 1990 and 1999, and investigated to which extent a 
temporal concentration trend can be observed after deregulation, by 
focusing on the mechanism of wave systems. In an hub and spoke 
traffic network, airlines typically operate synchronized waves of flights 
through these hubs. The aim of such a wave-system structure is to op- 
timize the number and quality of connections offered in the attempt 
to minimize waiting times for indirect connections. By comparing 
the flows of departing and arriving flights in an airport with an ideal 
type connection wave (given some transfer times), authors were able 
to identify the presence of waves and study how the structure and 
number of waves in an airport changed during the 1990s, when the 
air transport network became more and more similar to an hub and 
spoke system. Authors concluded that in the late 1990s major hubs 
had a clear wave system in place. For example, in Paris CDG air- 
port six clear daily waves could be identifled. The presence of these 
waves was then used to assess their effects on the quality of indirect 
connectivity. By introducing a suitable airport metric measuring the 
number of indirect connections weighted by transfer time and a rout- 
ing factor0, Burghouwt et al. [66J concluded that in 1999 only few 
airline hubs are highly competitive in the indirect market (essentially 
Frankfurt, Paris CDG, London Heathrow, and Amsterdam). More- 
over they showed evidence that European airlines have increasingly 
adopted wave-system structures or intensified the existing structure 
during the 1990s. Finally, they found evidence that, given a certain 
number of direct flights, airports adopting a wave-system structure of- 
fer generally better indirect connections than airports without a wave- 

^The routing factor is the ratio between the actual in-flight time indirect connection 
and estimated in-flight time direct connection based on great circle distance. It measures 
the quality of the indirect connection. 
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system structure. 

In order to assess the role of the European hub and spoke network, 
Mahghetti et al. | i67j investigated different notion of centrahty in the 
flight network, by considering the point of view of a passenger that 
must reach a destination in the shortest possible time. They first 
considered two purely topological measures. The first is the average 
shortest path length for an airport i, defined as the average of the 
minimum number of flights needed to reach an airport j from airport i. 
The second metric is the betweenness of an airport i, i.e. the number of 
minimal paths that passes through airport i. Since there are typically 
many shortest paths connecting two airports (and passing through 
z), authors also defined the essential betweenness as the number of 
unavoidable minimal paths passing through an airport i. An high value 
of the essential betweenness indicates that the airport is a bottleneck 
for the traffic in the system. Analyzing the European air transport 
network in 2007, they found a large heterogeneity of average shortest 
path length and different values for an airport, both if one considers 
only European flights and if one considers also world flights. This 
result is somewhat expected and typically differentiates hubs for main 
airlines from hubs for low cost carriers. 

Authors claim that, despite being interesting for the characteriza- 
tion of the centrality properties of the flight network, these topological 
measures do not help much in assessing the passengers' needs. In fact, 
a short path (in number of flights) can be useless for a passenger if 
the composing flights are unfrequent or if their scheduling makes im- 
possible to do the connection. For this reason Malighetti et al. [67] 
considered a speciflc day and for each pair of airports they calculated 
the shortest travel time between them for a passenger who wants to 
leave at a speciflc time and to arrive at destination within the same 
day. They used scheduled flight time data and allowed for at least 
one hour for a flight connection. Finally they considered the mini- 
mum among the possible starting time within the day and deflned the 
optimal path from the two airports as the one having this minimum 
travel time. This value is a passenger oriented metric for the distance 
between two airports in a given day. 

It is clear that this metric takes implicitly into account the geogra- 
phy of Europe. In fact, most of the airports with the smallest average 
optimal path time are also those in the central part of Europe, in part 
simply because flights are shorter from there than from a peripheral 
airport. Authors compared also the connectivity offered by three big 
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alliance networks (Oneworld, Sky Team, and Star Alliance) with that 
of the overall network. They found that roughly two thirds of the 
fastest indirect connections are not operated by the alliance system. 
The interpretation provided is that this could be exploited to enable 
a new passenger strategy of "self-help hubbing" . However, as authors 
warned, this result does not take into account properly travelers' util- 
ity, since the analysis is focused on time and neglects other important 
variables, such as prices, number of flights, loyalty programs, etc. 

A similar, but mathematically more sophisticated, approach to the 
indirect connectivity problem has been used in Ref. |68] where sched- 
uled networks have been introduced. A scheduled network for the 
air transport system is constructed starting from the ordinary flight 
network and by adding additional nodes on the links connecting two 
airports. The number of additional (secondary) nodes on a link is 
proportional to the traveling time needed to travel in the route asso- 
ciated to the link. The full network composed by primary (airports) 
and secondary nodes allows to compute the real time needed to go 
from an airport to another one, even if they are only indirected con- 
nected. The main advantage of this approach is that one can adapt 
several standard network metrics, such as mean path length, giant 
component, clustering coefficient, and tolerance to errors and attacks, 
to scheduled networks in order to take into account real travel time 
(and not just topological distance) as a metric and, more important, 
to take properly into account indirect connectivity. Authors then ap- 
plied their method to a small sample of flight data of 40 European 
airports and measured the efficiency of the network as a function of 
the time of the day. As expected, they found a very low efficiency 
during the night, meaning that passengers should wait in the airport 
for some connections until the next morning. Moreover, they identi- 
fied three peaks of connectivity during daytime, corresponding with 
the moments of high network connectivity (see Figure 5). 

4.2 Air traffic jams 

The commercial airline trafhc is made of scheduled flights. However 
due to several possible reasons, such as adverse weather conditions, 
operational problems, and high traffic volume, the actual dynamics 
of flights can be drastically different from the scheduled one. Due to 
the strong interconnectedness of the air traffic system, it is likely that 
deviations from the scheduling propagate in space and time, i.e. in 
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Figure 5: Evolution of three metrics for the 40 busiest European airports, at 
different hours of the day. The metrics are the proportion of active connec- 
tions (black solid line, left scale), the efficiency of the network (blue dashed 
line, left scale), and the minimum rotation time, i.e., the minimum time re- 
quired for a crew to return to the departure airport (red dotted line, right 
scale). Reprinted with permission from Ref. [68]. ©2009 by the American 
Institute of Physics. 



other airports or routes in the near future. A big engineering challenge 
is to design the whole system in such a way to be resiUent to these 
shocks, i.e. to be able to return quickly to a normal state after a shock. 
Moreover air traffic system is a specific instance of traffic where jams 
can appear for no apparent reasons. A small set of recent papers 
started to investigate theoretically and empirically this problem. 

In Ref. [69j, for example, authors proposed a network-based model 
of the air transport system that simulates the effect of traffic dynam- 
ics and shows the appearance of jams. Specifically, in the model a 
random (Erdos-Renyi) network describes the topology of a set of in- 
terconnected airports. Each airport is characterized by an exogenously 
given capacity, i.e. the maximal number of aircraft per unit time that 
the airport can handle in an ideal situation. This ideal capacity is per- 
turbed (diminished) by a random noise term. Moreover each link of 
the network is weighted and the weight measures the number of time 
steps that an aircraft needs in order to complete the route. A series 
of simple queuing rules describes the interplay between incoming and 
out-coming aircraft flows in an airport. Simplifying a bit, if the input 
flow is larger than the capacity, the output flow will be equal to the 
capacity and the other aircraft will remain in a queue. Only when 
the input flow becomes smaller than the capacity it will be possible to 
remove the waiting aircraft from the queue. The model is simulated 
with Monte Carlo methods. 
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Figure 6: Phase diagram relating the percentage of aircraft not stuck in a 
node's queue P as a function of the network aircraft density, for the European 
air transport network composed of 858 nodes. Reprinted with permission 
from Ref. |69]. ©2009 by Elsevier. 




The key system indicator, P, is the percentage of aircraft that 
are not stuck in a node's queue, measured in the steady state. Note 
that P actuahy measures the network's efficiency as far as it gives 
the flow rate which is difl^using as compared to the flow rate which is 
stuck. The key flnding of Ref. [69] is that by increasing the aircraft 
density (number of aircraft), the system undergoes a phase transition. 
This is testifled by the fact that the expected value of P sharply 
deviates from the efficient phase P = 1 when the aircraft density is 
larger than a threshold. After this threshold P declines, as expected. 
Correspondingly, the variance of P, goes abruptly from zero to an high 
value when the aircraft density is larger than the threshold described 
above. These and other evidences indicate that the system undergoes a 
jamming transition, similarly to what observed in other traffic systems 
|70j . One may ask whether network topology plays a role in this 
observation. Authors considered the topology of the (real) European 
air traffic system composed by 858 airports and 11170 flights and they 
found qualitatively the same result, i.e. the emergence of a jamming 
phase transition for a given value of the aircraft density (see Figure 




6). 
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4.3 Epidemic spreading 



Another recent stream of research has considered the role of air trans- 
portation network in the propagation of global epidemics and in the 
assessment of its predictability. The majority of "classical" models of 
epidemic spreading considers a set of individuals located very closely 
one to each other, so that the connections between individuals that 
can propagate the epidemic are short range (in space). However air 
transport network provides a mechanism for long-range heterogeneous 
connections that can change dramatically the diffusion properties of 
an epidemics. 

In Ref. |71] Colizza et al. developed a model of a set of more than 
3000 large cities worldwide where a major airport is present. They 
use lATA data on the number of available seats on any given flight 
connection for the year 2002 and they complemented this dataset with 
census data on the population of the large metropolitan area served 
by each airport. Interestingly, they found that the number of pas- 
sengers (seats) scales as the square of the population. For each city, 
they simulated a classical SIR epidemic model, where each individual 
is either susceptible (S), or infected (I), or recovered (R). Then they 
used the air traffic flow data to simulate the mobility of individuals 
from an airport to another. In this way they were able to simulate 
the epidemic spreading in the world and to assess the role of the air 
transportation network for the global pattern of emerging diseases. 
To this end they compared the simulations calibrated on the real air 
traffic data with simulations for two benchmarks. In the first one, the 
air traffic network is a random Erdos-Renyi graph, while in the second 
one the network is a graph with the same topology of the real system 
but fluxes and populations are taken as uniform and equal to the av- 
erage value of the corresponding real variables. In this way authors 
could assess the relative importance of topological (network) and met- 
ric (population, SIR rates, etc) variables by comparing simulations of 
the two benchmark cases to simulations of the model fully calibrated 
on real data. The striking result is that the model where the topology 
of the air network is faithfully reproduced but the weights are homo- 
geneous shows simulations much more similar to those of the model 
fully calibrated on real data than the simulations of the random graph 
model. This finding strongly indicates that (at least in the framework 
of the model) "the air-transportation-network properties are responsi- 
ble for the global pattern of emerging diseases" [7T]. The same result 
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holds if one consider the predictabihty of epidemic spreading. By 
defining a suitable measure of the sensitivity of simulations to ran- 
dom perturbation, Colizza et al. [7l] found that the predictability (or 
more precisely, the reproducibility or the model's sensitivity) of the 
random graph model is much higher than the model fully calibrated 
on real data. Again, the model where only the real air transport net- 
work topology is reproduced displays patterns similar to the model 
fully calibrated on real data. These results point out that topology 
of the air transport network have an important role, not only for the 
mobility of people, but also for the dynamics of entities that depend 
on human mobility. 

5 Resilience and vulnerability 

In this last Section we review some results on the resilience of the air 
transport network, i.e., its ability to adjust its functioning prior to, 
during, and following internal and external disturbances, so that it can 
sustain required operations under both expected and unexpected con- 
ditions |72]. For example, few preliminary studies (for example, |51] ) 
found a positive correlation between node's topological properties and 
typical fraction of delayed flights. In spite of its relevance for passen- 
gers and society in general [17J however, little effort has been devoted 
to the understanding of the relationships between the topology of the 
air transport networks and the vulnerability of its dynamics. 

Two significant exceptions can be found in the literature. The first 
one, by Chi and coworkers [73], analyzes how the main topological 
properties of the US air transport network are changed by random 
failures and attacks. The former effect is analyzed by de-activating 
airports at random, and thus simulating random disturbances like 
emergency situations or adverse weather. The latter is investigated 
by deleting the most connected nodes (as in an intentional terrorist 
attack). As known in more general settings |74j . scale- free networks 
are extremely resilient to random failures. However this comes at 
a high price, because they are also extremely vulnerable to targeted 
attacks. The de-activation of the 10% of the most connected airports 
is enough to reduce the topological efficiency [751 of the network of 
25%. 

The resilience of the air transport system to random failures, as 
obtained in Ref. [73j, is partly aligned with the experience that we 
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all have as passengers. In spite of multiple failures that may appear 
in small airports across the network, the dynamics of the system as 
a whole is seldom disturbed. Yet, a complementary problem is rep- 
resented by those events that push the dynamics of the system far 
away from its normal point of operation. For instance, black swans as 
large strikes or the eruption of a volcano dramatically affect the per- 
formance of the system. The eruption of the Eyjafjallajokull volcano 
in 2010 is a clear example of such random events that have larger than 
expected consequences. 

This problem has been tackled in Ref. [7B]. Specifically, a set of 
random networks has been created, in which the main properties of the 
topology (i.e., its scale- free nature) and the spatial position of nodes 
have been maintained as found in the European air transport network. 
Results indicate that the severe disruptions observed in 2010 are ex- 
plained by the geographical correlation of the disturbances (which are 
not, thus, completely random, as in the model of Ref. ^3]), and by 
the geographical correlation of hubs, which concentrates in the centre 
of Europe. The proposed solution is to move some hub airports from 
Germany to peripheral regions. However, although this may improve 
the resilience against black swans, the economical consequences for 
airlines would probably exceed the expected benefits. 

6 Conclusions and open lines of re- 
search 

In conclusions, we have presented a short review of the recent use of 
complex network methods for the characterization of the structure of 
air transport and of its dynamics. We have shown that most of the 
published researches have focused on the topological and metric prop- 
erties of flight networks, where nodes represent airports, and links are 
created between pairs of them, bringing information on the presence 
and the frequency of flights. Speciflcally, these papers can be classi- 
fied within three main families. Firstly, some of them propose a simple 
characterization of the topology of the networks, without considering 
their evolution through time. The recent change from a point to point 
to an hub-and-spoke system has triggered a series of studies, aimed at 
identifying and characterizing this transition, by monitoring the evo- 
lution of network's characteristics through time. Finally, some works 
have focused on the dynamics on the network, as, for instance, on the 



23 



movement of passengers and the epidemic spreading. However, we be- 
lieve that this is just a starting point of a fruitful breeding between air 
transport science and complex network theory. Many different types 
of networks can be defined by taking into account variables or phe- 
nomena that have not been investigated with networks. Therefore, we 
anticipate that many interesting contributions will be published in the 
near future about air transport networks involving, for example, air- 
ways and navpoints, delays, safety events, crews and physical aircraft, 
sectors, etc. 

In the near future our society is likely to face a significant growth 
of air traffic. Radical organizational and institutional changes are al- 
ready taking place to accommodate this increase, for instance with the 
progressive integration of the still nationally fragmented Eu airspace 
management, thanks to the Single European Sky initiative and the 
SESAR program. The current system will drastically change, thus 
making a primary research area of the development of innovative net- 
work management methods and tools. We believe that complex net- 
work theory could give a significant contribution to this challenge. 

7 acknowledgement 

Authors acknowledge S. Micciche, R. N. Mantegna, and S. Pozzi for 
useful discussion. This work has been developed as part of the activ- 
ities of the Complex World Network (www. complcxworld.eu). Work 
presented therein was co-financed by EUROCONTROL on behalf of 
the SESAR Joint Undertaking in the context of SESAR Work Package 
E, project ELSA. The paper reflects only the authors' views. EURO- 
CONTROL is not liable for any use that may be made of the infor- 
mation contained therein. 

References 

[1] Air Transport Action Group, The economic and social benefits of 
air transport 2008, (2008) 

[2] T. Bolic, Z. Sivcev, Transportation Research Record 11-3133, 

(2011) 136-143 

[3] M. Mazzocchi, F. Hansstein, M. Ragona, CESifo Forum 11, (2010) 
92-100 



24 



[4] S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, D.-U. Hwang, 
Physics Reports 424, (2006) 175-308 

[5] R. Albert, A.-L. Barabasi, Rev. Mod. Phys. 74, (2002) 47-97 

[6] L. da F. Costa, F. A. Rodrigues, G. Travieso, P. R. Villas Boas, 
Advances in Physics 56, (2007) 167-242 

[7] L. da F. Costa, O. N. Oliveira, G. Travieso, F. A. Rodrigues, P. 
R. V. Boas, L. Antiqueira, M. P. Viana, L. E. C. Rocha, Advances 
in Physics 60, (2011) 329-412 

[8] F. Liljeros, R. Edhng, L. A. N. Amaral, H. E. Stanley, Y. Aberg, 
Nature 411, (2001) 907-908 

[9] R. Pastor-Satorras, A. Vazquez, A. Vespignani, Physical Review 
Letters 87, (2001) 258701 

[10] E. Bullmore, O. Sporns, Nature Reviews Neuroscience 10, (2009) 
186-198 

[11] O. Sporns, D. R. Chialvo, M. Kaiserc, C. C. Hilgetag, Trends in 
Cognitive Sciences 8, (2004) 418-425 

[12] P. Crucitti, V. Latora, S. Porta, Physical Review E 73, (2006) 
036125 

[13] S. Porta, P. Crucitti, V. Latora, Physica A 369, (2006) 853-866 

[14] P. Sen, S. Dasgupta, A. Chatterjee, P. A. Sreeram, G. Mukherjee, 
S. S. Manna, ar Xiv:cond-mat /0208535 | y2 

[15] V. Latora, M. Marchiori, Physica A 314, (2002) 109-113 

[16] P. Angeloudis, D. Fisk, Physica A 367, (2006) 553-558 

[17] European Commission, Roadmap to a Single European Trans- 
port Area - Towards a competitive and resource efficient transport 
system, (2011) 

[18] M. Kurant, P. Thiran, Physical Review Letter 96, (2006) 138701 

[19] L. Adamic, Lecture Notes in Computer Science 1696, (1999) 443- 
452 

[20] R. Albert, H. Jeong, A.-L. Barabasi, Nature 401, (2001) 130-131 
[21] A. Vazquez, R. Pastor-Satorras, A. Vespignani, 



arXiv:cond-m at /0206084vl 
[22] D. D. Han, J. H. Qian, J. G. Liu, |arXiv:physics/0703193| y2 
[23] W. Li, X. Cai, Physical Review E 69, (2004) 046106 



25 



[24] Y.-Q. Qu, X.-L. Xu, S. Guan, K.-J. Li, S.-J. Pan, C.-G. Gu, 
Y.-M. Jiang, D.-R. He, larXiv:1010.2572l (2010) 

[25] X.-L. Xu, Y.-Q. Qu, S. Guan, Y.-M. Jiang, D.-R. He, Europhysics 
Letters 93, (2011) 68002 

[26] K.-Q. Cai, J. Zhang, W.-B. Du, X.-B. Cao, Chin. Phys. B 21, 
(2012) 028903 

[27] N. Pyrgiotis, K. M. Malone, A. Odoni, Transportation Research 
Part C, in press 

[28] N. Bania, P. W. Bauer, T. J. Zlatoper, Transportation Research 
Part E 34, (1998) 53-74 

[29] F. Lillo, S. Pozzi, A. Tedeschi G. Ferrara, G. Matrella, F. Lieu- 
taud, B. Lucat, A. Lieu, Coupling and Complexity of Interaction of 
STCA Networks. Proceedings on the EUROCONTROL Conference 
8th Innovative Research Workshop & Exhibition, Bretigny-sur-Orge 
(France) December 1-3 2009. 

[30] M. Alderighi, A. Cento, P. Nijkamp, P. Rietveld, Transport Re- 
views 27, (2007) 529-549 

[31] Y. H. Chou, Transportation Planning and Technology 14, (1990) 
243-258 

[32] M. E. O'Kelly, H. J. Miller, Journal of Transport Geography 2, 
(1994) 31-40 

[33] S. Berry, M. Carnall, P. T. Spiller, NBER Working Paper No. 
5561, (1996) 

[34] R. Guimera, S. Mossa, A. Turtschi, L. A. N. Amaral, PNAS 102, 
(2005) 7794-7799 

[35] A. Barrat, M. Barthelemy, R. Pastor-Satorras, A. Vespignani, 
PNAS 101, (2004) 3747-3752 

[36] A. Barrat, M. Barthelemy, A. Vespignani, J. Stat. Mech. 2005, 
(2005) 

[37] C. Li-Ping, W. Ru, S. Hang, X. Xin-Ping, Z. Jin-Song, L. Wei, 
C. Xu, Chin. Phys. Lett. 20, (2003) 1393-1396 

[38] Z. Xu, R. Harriss, GeoJournal 73, (2008) 87-102. 

[39] J. Wang, H. Mo, F. Wang, F. Jin, Journal of Transport Geogra- 
phy 19, (2011) 712-721 

[40] G. Bagler, Physica A 387, (2008) 2972-2980 



26 



[41] M. Sapre, N. Parekh, Lecture Notes in Computer Science 6744, 
(2011) 376-381 

[42] M. Guida, F. Maria, Chaos Solitons and Fractals 31, (2007) 527- 
536 

[43] J. Quartieri, M. Guida, C. Guarnaccia, S. DAmbrosio, D. 
Guadagnuolo, International Journal of Mathematical Models and 
Methods in Applied Sciences 2, (2008) 312-316 

[44] M. Zanin, J. M. Buldu, P. Cano, S. Boccaletti, Chaos 18, (2008) 
23103 

[45] A.-L. Barabasi, Science 24, (2009) 412-413 

[46] L. C. Freeman, Sociometry 40, (1977) 35-41 

[47] P. Erdos, A. Renyi, Publications of the Mathematical Institute 
of the Hungarian Academy of Sciences 5, (1960) 17-61 

[48] M. P. H. Stumpf, C. Wiuf, R. M. May, PNAS 102, (2005) 4221- 
4224 

[49] Z. Wu, L. A. Braunstein, V. Colizza, R. Cohen, S. Havlin, H. E. 
Stanley, Physical Review E 74, (2006) 056104 

[50] D. D. Han, J. H. Qian, J. G. Liu, Physica A 388, (2009) 71-81 

[51] F. Lillo, S. Micciche, R.N. Mantegna, V. Beato, S. Pozzi, ELSA 
Project: Toward a complex network approach to ATM delays anal- 
ysis, Proceedings on the INO 2011 Conference, Toulouse (France, 
2011) 

[52] M. E. J. Newman, Physical Review E 70, (2004) 056131 

[53] S. D. Barrett, Transportation 16, (1990) 311-327. 

[54] J. Berechman, J. de Wit, Journal of Transport Economics and 
Politics 3, (1996) 251-274 

[55] O. Betancor, J. Campos, Proceedings of the 1997 Air Transport 
Research Group (ATRG) of the WCTR Society, (1997). 

[56] M. B. Rey, Journal of Air Transport Management 9, (2003) 195- 
200 

[57] G. Burghouwt, J. Hakfoort, Journal of Air Transport Manage- 
ment 7, (2001) 311-318 

[58] F. Jin, F. Wang, Y. Liu, The Professional Geographer 56, (2004) 
471-487 



27 



[59] J. Wang, F. Jin, Eurasian Geography and Economics 48, (2007) 
469-480 

[60] X. Ma, M. F. Timberlake, GeoJournal 71, (2008) 19-35 

[61] J. Zhang, X.-B. Cao, W.-B. Du, K.-Q. Cai, Physica A 389, (2010) 
3922-3931. 

[62] L. E. C. Rocha, J. Stat. Mech., (2009) P04020 

[63] F. Dobruszkes, Journal of Transport Geography 14, (2006) 249- 
264 

[64] A.R. Goetz and C. J. Sutton, Journal of the Association of Amer- 
ican Geographers 87 (1997) 238-263. 

[65] K. Button, Journal of Air Transport Management 8 (2002) 177- 
188. 

[66] G. Burghouwt and J. de Wit, Journal of Air Transport Manage- 
ment 11 (2005) 185-198. 

[67] P. Malighetti, S. Paleari, and R. Redondi, Journal of Air Trans- 
port Managment 14, (2008) 53-65. 

[68] M. Zanin, L. Lacasa, M. Cea, Chaos 19, (2009) 023111 

[69] L. Lacasa, M. Cea, and M. Zanin, Physica A 388 (2009) 3948- 
3954. 

[70] D. Helbing, Review of Modern Physics 73 (2001) 1067-1141. 

[71] V. Colizza, A. Barrat, M. Barthelemy, and A. Vespignani, Pro- 
ceedings of the National Academy of Sciences USA 103 (2006) 
2015-2020. 

[72] E. HoUnagel, J. Paris, D. D. Woods and J. Wreathall, Resilience 

Engineering in Practice (Ashgatc, 2011). 

[73] L. P. Chi, X. Cai, International Journal of Modern Physics B 18, 
(2004) 2394 2400. 

[74] R. Albert, H. Jeong, A.-L. Barabasi, Nature 406, (2000), 378-382 

[75] V. Latora, M. Marchiori, Europ. Phys. Journ. B 32, (2003) 249- 
263. 

[76] S. M. Wilkinson, S. Dunn, S. Ma, Nat Hazards 60, (2012) 1027- 
1036. 



28 



